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Abstract 

Hadron production in lepton-nucleus deep inelastic scattering is studied in an ab- 
sorption model. In the proposed model, the early stage of hadronization in the 
nuclear medium is dominated by prehadron formation and absorption, controlled 
by flavor-dependent formation lengths and absorption cross sections. Computations 
for hadron multiplicity ratios are presented and compared with the HERMES ex- 
perimental data for pions, kaons, protons and antiprotons. The mass-number depen- 
dence of hadron attenuation is shown to be sensitive to the underlying hadronization 
dynamics. Contrary to common expectations for absorption models, a leading term 
proportional to A 2 / 3 is found. Deviations from the leading behavior arise at large 
mass-numbers and large hadron fractional momenta. 
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1 Introduction 



Recent HERMES results give precise data on hadron production in deep in- 
elastic scattering (DIS) of 27.6 GeV positrons on D, He, N, Ne, and Kr nuclei 
[1-3]. The main observable is the multiplicity ratio, Rm, defined as the ratio 
between the hadron multiplicity on nucleus and on deuterium. Rm has been 
studied as a function of the hadron fractional momentum z, of the virtual 
photon energy u, of its virtuality Q 2 , for different hadrons and for different 
nuclei. The use of nuclear targets allows to study the hadronization process 
near the interaction point of the photon and probe hadron formation only 
some fermi away. 

There are many proposed models to describe how hadronization in presence 
of a nucleus evolves in space and time. The model computations range from a 
color string breaking mechanism [4-9] with final state interaction [10, 11], to 
gluon bremsstrahlung for leading hadrons [12, 13] and pure energy loss mod- 
els [14-16]. Two classes of models of hadron formation compete with each 
other. The first one is based on nuclear absorption, where the color of the 
struck quark is neutralized after a short time by the formation of a prehadron, 
the predecessor of the final hadron, which interacts with surrounding nucleons 
and is absorbed on its way out of the nucleus. The second model assumes 
energy loss from medium-induced radiation of the struck quark as dominant 
process. The first interpretation emphasizes the hadronic aspects of particle 
production, the second one focuses on partonic degrees of freedom and post- 
pones hadronization outside the nucleus. Neglecting the production of sec- 
ondary particles, both mechanisms reduce the number of hadrons emerging 
from the nucleus. The current belief is that the dependence of hadron atten- 
uation from the mass-number A of the target nucleus can differentiate the 
two processes: in the absorption model, hadron attenuation is commonly be- 
lieved to be proportional to the path length of the (pre)hadron in the nucleus 
(oc A 1 / 3 ), whereas in the energy loss model the attenuation is supposed to de- 
pend on the square of the distance the quark traverses in the nucleus (oc A 2 / 3 ). 
A careful measurement of the A-dependence of the nuclear attenuation would 
therefore allow to discriminate between the two different processes. 

In a previous paper [17], we have calculated the nuclear modifications of 
hadron production in DIS in an absorption model. The hadron formation 
length has been computed analytically in the framework of the Lund string 
model [5, 18] as a two-step process. In the first step a quark-antiquark pair 
from the break-up of the color string forms a prehadron. In the second step 
the final state hadron is created. Both the prehadron and the hadron, if they 
form inside the nucleus, interact with target nucleons and may be absorbed. 
Setting the prehadron nucleon cross section equal to the hadron nucleon cross 
section in Ref. [17] we had to increase the formation length by using an ef- 
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fective string tension k—OA GeV/fm, much smaller than the expected k ~ 1 
GeV/fm (see Ref. [19] for example). 

In this paper we correct the model by combining a realistic formation length 
based on the expected string constant with smaller prehadronic cross sections. 
We pay special attention to the flavor dependence of the formation length, 
which is naturally induced by the Lund model, and calculate the multiplicity 
ratio R M for different hadron species h as functions of kinematic variables z 
and v. Finally, we study the A-dependence of the nuclear attenuation 1 — R M , 
and challenge common expectations by showing, analytically and numerically, 
a A 2 ^ dependence for the absorption model. 

The paper is organized as follows. In section 2, the new developments of the 
model are described. Section 3 is devoted to the analytic calculation of the 
A-dependence of hadron attenuation. In Section 4, numerical computations 
of R M and of the A-dependence of 1 — Rm are presented and compared with 
experimental data. Section 5 presents a summary and our conclusions. 



2 Development of the model 



In this section we report the main developments of the model, which has been 
extensively described in Ref. [17]. The experimental multiplicity ratio R h M 
between nucleus and deuterium is a double ratio and can be measured as a 
function of z, v or Q 2 . As a function of z it is given by: 



* ^ . - dN\{z) I 1 dN h D {z) 



R M(Z) - N DIS dz I N DIS dz ■ 

The upper ratio is the number of produced hadrons of species h with energy 
fraction z, normalized to the total number of deep inelastic scattering events 
on a nuclear target with mass number A. The lower ratio is the corresponding 
expression for a deuterium target. The theoretical calculation of the hadron 
multiplicity on a nucleus A is based on the computation of the hadron survival 
probability Sf h and the fragmentation function D'j in the nucleus: 



1 dN\(z) 
N® IS dz 



its 

da lq 



dx dv 
exp. cuts 



^J2e}qf(x,UQ 2 )^Sf A (z,u)D h f (z,UQ 2 ). (2) 

The total lepton-nucleus DIS cross section a lA is calculated in leading order by 
an integration over the parton distribution functions g/ and the lepton-quark 
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cross section da lq /dxdv, including the experimental cuts: 



a lA = j dxdu^e) q f{ X ,UQ 2 )^- ■ (3) 
exp. cuts 



The theoretical calculation in Eqs. (2-3) takes into account several nuclear 
effects. The partial deconfinement model [20] is used to express the nuclear 
parton distribution function and the nuclear fragmentation function: both 
the free-nucleon parton distribution function qf(x,Q 2 ) and the free-nucleon 
fragmentation function Dj(z,Q 2 ) are rescaled with a factor £4 = £a(Q 2 ) due 
to the hypothesis that quarks in a bound nucleon have access to a larger region 
in space than in free nucleons, i.e. 

A a > Ao (4) 

where A^ is the confinement scale of a bound nucleon, which is assumed to be 
proportional to the overlap of nucleons inside the given nucleus, and A is the 
confinement scale of a free nucleon. Effectively, the DGLAP evolution of the 
nuclear structure functions and the nuclear fragmentation functions covers a 
larger interval in momentum compared with the corresponding functions in 
the nucleon at the same scale Q 2 which implies an increased gluon shower. 
Details are presented in [17]. 

The hadron survival probability S^ h represents the probability that the hadron 
emerges from the nucleus without having interacted with the nucleus during 
its evolution from the quark to its final state. The survival probability depends 
on the struck quark flavor /, the hadron species h and the target nucleus A. 
For simplicity, any scattering of the prehadron with a target nucleon is as- 
sumed to lead to absorption. In principle the observed hadron with energy 
fraction z might have undergone several inelastic rescatterings before being 
absorbed, whereby its z is degraded. However, these multiple scattering pro- 
cesses are suppressed at large z since the fragmentation function is rapidly 
falling at z — > 1, and are expected to be negligible at z > 0.3 — 0.4. In order to 
calculate Sf h we modify the simplistic Bialas-Chmaj [6] absorption formulae 
used in Ref. [17], which considers the hadronization process to be a decay 
process, by fully coupling the evolution of the quark into a prehadron and a 
hadron with the absorption processes. The quark decays into the prehadron 
and has an average lifetime which equals the average formation length (/*) 
of the prehadron after the interaction of the virtual photon 7* with a quark 
q. The prehadron itself decays into the hadron and has an average lifetime 
(A/) given by the difference of the average hadron and the average prehadron 
formation lengths: 

(A/) = (i h ) - <r> (5) 

Here we omit for ease of notation the dependence of the prehadron and hadron 
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formation length on the struck quark flavor /, the hadron species h, the energy 
fraction z and the virtual photon energy v. We will come back later on this 
topic. We compute the formation lengths in the framework of the standard 
Lund string fragmentation model. The Lund model is a semiclassical model 
which provides a formation length distribution of the produced prehadrons 
and hadrons derived from classical relations among the production points of 
the hadrons without quantum fluctuations. Since, on the other hand, string 
breaking is a quantum process akin to quantum tunneling, we can expect the 
Lund model to be a good ansatz for the full probability distribution only up 
to the first few moments. For the sake of simplicty, in this paper we consider 
only the first moment, i.e., the average formation length. 

If the initial 7*g interaction occurs at longitudinal coordinate y, the probabil- 
ities that the intermediate state at time y' is a quark, P q (y,y'), a prehadron, 
P*(y, y')i or a hadron, P h (y,y'), satisfy the following differential equations: 

dP q (y,y') _ P q {y,y') 

By' (I*) 
dP*(y,y') _P q (y,y') P*(y,y') P*(y,y') 



dy> ~ (I*) 
dP h (y,y') _ P*(y,y') 
dy> ~ (A/) 



with initial conditions 



P q (y,y' = y)=l 

P*(y,y' = y) = 

P h (y,y' = y) = 0. 



(A/) A*(y') 



(6) 



The mean free path of the prehadron K(y') and the mean free path of the 
hadron \h(y') are: 

K h (y') = An (7) 

Ap A (y')cr*,h 

where pa is the nuclear density normalized to unity and h are the respective 
cross sections of a prehadron and of a hadron. The dependence of pa on the 
impact parameter b is suppressed to simplify the notation. Motivated by the 
experimental fits in Ref. [1,21] we neglect any final state interaction of the 
struck quark with the nuclear environment (cf. Ref. [17]). The solutions for 
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the probabilities P q ,P*, Ph can be obtained analytically: 

P q (y,y')=e-te (8) 



y 



P*(y,y') = y dx— -e W e - (9) 



V X a; — y p a: —a: C j A I \ 

f r e ~ a * J dsA PA{s) e "(siy -°"h J dsAp A (s) 
P h (y,y') = j dx' j dx-^-e - "(A/)" 6 ' ^ 
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There is a connection between the resulting probability distribution P q (y,y') 
to find a quark at position y' if the initial interaction took place at y and the 
Lund model prehadron formation length distribution. The prehadron forma- 
tion length distribution in our model is simply given by the first derivative of 
P q with respect to y' with an opposite sign. Its first moment equals the first 
moment of the Lund distribution. Its second moment deviates from the Lund 
distribution's second moment by less than 20% for z > 0.5, with larger devia- 
tions at smaller z. An analogous analysis also holds for the hadron formation 
length distribution. We define the survival probability Sf h of the hadron as 
the probability that the hadron is not absorbed by the nucleus. It is obtained 
as linij/^oo Ph(y, y') after integration over all interaction points and impact 
parameters: 



OO X 



Sf h = J d 2 b J dypA(b,y) J dx' J dx 



-oo y y 



e JFJ -<J* J dsAp A (s) g JaiJ —<*h J dsAp A (s) 

X /, ; e x , A e *=' . (11) 

(I*) (A/) v 1 



The formation lengths are computed in the framework of the Lund model [5, 
18], for which we show a schematic space-time picture in Fig. 1. Due to the 
relatively high measured values of u, hadronization is supposed to occur near 
the light cone which justifies the neglection of all masses. This is definitely 
true for hadrons with large z. The target nucleons Ni,N 2 are at rest in the 
laboratory frame. Since the struck quark and the leading hadrons move close 
to the light cone, we have enlarged in our illustration the negative light cone 
momenta and diminished the positive ones for clarity. For the same reason 
we have omitted a part of the space-time evolution which is indicated by the 
two dashed vertical lines on the left of the figure. The hadrons Hi are formed 
at the points F$ and the index % gives their rank, counted from right to left. 
Each point Pj denotes the production point of one quark anti-quark pair due 
to string breaking. We exemplarily consider the production point P 2 in the 
following. Since the produced antiquark is still connected via a string to the 
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Fig. 1. Schematic space time picture of hadronization in the Lund model. 



struck quark q , we associate with the combined object (qqo) the prehadron 
if 2 evolving into the hadron H2. Therefore, the "production" points Pi are the 
formation points of the i-th rank prehadron. The distance L from the struck 
quark turning point to the j*q interaction point sets the scale of hadronization 
and is proportional to the energy transferred to the quark: 



v 

K 



;i2) 



For the calculation in this paper we use the expected vacuum string constant 
[4,7,9,19] 



k = 1 GeV/fm . 



(13) 



In order to be consistent with the partial deconfinement model, we rescale 
the vacuum string constant in the nucleus. Since the string constant is the 
physical quantity that sets the confinement scale, a larger confinement scale 
corresponds to a smaller string tension in the nucleus: 



Ai 



(14) 
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In addition, we explicitly consider the dependence of the formation time on 
hadron species. The rank 1 hadron has the struck quark as one of its con- 
stituents. Therefore, the type of particles which can be produced as a rank 1 
depends on the flavor of the struck quark. As an example, in Fig. 2 we show 
the different processes which lead to the production of a K + and K~ meson 
from a struck u quark. The positive kaon can be directly formed at rank 1, 
whereas K~ needs a least a second quark pair produced, i.e. it can only be 
produced from rank 2 on. More in general, a negative kaon (K~ = su), as 
well as an antiproton {p = uud), cannot be formed as a first rank particle 
by a struck valence quark in the nucleon. As a first rank particle, the kaon 
and the antiproton can only be formed from a struck sea quark, which is a 
subdominant process at HERMES (q u ,d/c[s ~ 0.05) included in our calcula- 
tion. In any case, K~ and p can be formed from quarks inside the color string 
,i.e. as higher rank particles. Therefore, kaons and antiprotons are dominantly 
produced as higher rank particles with shorter formation lengths. Hence, the 
Lund string model naturally produces a flavor-dependent formation time. In 
order to get analytic formulae for the average formation lengths eqns. (15) and 
(16) we make several approximations: we compute the average formation time 
with the standard Lund string fragmentation function f(u) oc (1 — u) Da (cf. 
Ref. [17]) which does not depend on the mass of the produced hadrons. We do 
not take into account the hadron masses in the determination of the relation 
between the hadron energy and the hadron production point. The main effect 
of an inclusion of hadron masses in the calculation is that the average forma- 
tion length is cut off at small values of z due to energy conservation but that 
the general behavior at larger z is not significantly affected. We average over 
an infinite number of produced hadron ranks. Furhermore, we do not take into 
account gluonic excitations of the strings and neglect transverse momenta in 
the computation of the average formation time. For details about the deriva- 
tion of the average formation times we refer to Ref. [17]. Here we only show 
how the mean formation time of a hadron which can be produced from all 
ranks n > 1 differs from the formation time of a hadron which can only be 
produced from rank 2 on, i.e. n > 2. The average prehadron formation length 
for hadrons which are producible as first rank particle reads: 



<*>!> 



l+Da 

l + C + (D a - C)z 



(l-z)zL 



l + ^^^^ 2 F l (2 + D a ,2 + D a ;3 + D a , Z - 



2 + D a z 2+D « 



(15) 



and the average prehadron formation length for particles which are not pro- 
ducible as a first rank particle is given by: 
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Fig. 2. Schematic picture describing the fragmentation of an up quark into K + and 
K~ . A similar picture works also for fragmentation into p and p, respectively. 
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Fig. 3. Computed prehadron formation lengths when an up quark is struck by the 
virtual photon. Left: When a 7r + , K + or p is observed, the corresponding prehadron 
can be created at rank n > 1. Right: When a tt~ , K~ or p is observed, the corre- 
sponding prehadron can be created only at rank n > 2. 



1-2)21 
1 + D a 1 



2 + £>„ z 2 + D ° 



2 F 1 (2 + D a ,2 + D a ;3 + D a ; 



(16) 



Here 2-F1 is the Gauss hypergeometric function, the parameters C and D a with 
a = (q, qq) arise from the string fragmentation function. Furthermore, these 
parameters select whether a meson (a = q) or a baryon (a = qq) is produced 
(see Ref. [17]). Their numerical values are given in Ref. [18] by C = D q = 0.3 
and D qq = 1.3. There is a simple rank-independent connection between the 
prehadron and the hadron formation length for fixed energy fraction z carried 
by the hadron : 



(*>w> = < Z >i,2> + *L ■ 



(17) 



The present model is based on string fragmentation in (l+l)-dimensions. Re- 
alistic hadrons also have a transverse extension. We expect that the prehadron 
does not yet have the full size of the hadron and therefore interacts with a 
smaller cross section, a* < Oh- We fit the prehadronic cross section to repro- 
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z z 

Fig. 4. Multiplicity ratios for K~ and p on krypton are shown as function of z in the 
left and right panel respectively. Solid (dashed) curves are computed with (without) 
hadron size rescaling in the formation length. The dotted curve shows the published 
result [17] for K~ without the above described model modifications. The data are 
taken from Ref. [2]. 



duce the pion multiplicity ratio on 84 Kr [2] and we obtain: 
2 

cr* = -<Th ■ ( 18 ) 

Such a prehadronic cross-section is in agreement with Ref. [11], where the 
multiplicity ratio of charged hadrons for the HERMES experiment is computed 
with a prehadronic cross-section increasing quadratically during the hadron 
formation time tj = 0.5 fm/c from 0.3cr/j to the full hadronic cross-section. 
On the average one obtains a prehadronic cross-section (a*) = 0.63(7^. We use 
the proportionality factor 2/3 for all other hadronic species. Empirically we 
know that baryons have a larger radius than 7r mesons, which are larger than 
K mesons. We renormalize the string tension, i.e. the quantity which sets the 
confinement scale to the confinement scales r\ of real hadrons taken from Ref. 
[22] 

2 

Kfc = K^f. (19) 
r h 

This produces an increased formation length for baryons due to the larger 
size of the proton compared with pions. A comparison between the computed 
formation lengths of pions, kaons, protons and antiprotons is shown in Fig. 3. 
The flavor dependence due to the mechanism exemplified in Fig. 2 decreases 
the average formation length of the negative mesons by about 30% compared 
to the positive ones, and by more than 40% for p compared to p. The effect 
of the baryon's larger size results in a slightly longer baryon formation time 
compared with the meson ones. 

The effect of the modified string tension Kh on the K~ and p multiplicity ra- 
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tios for krypton is shown in the left and right parts of Fig. 4, respectively. The 
dashed line represents the computation with the unrescaled string tension k, 
while the solid line shows the effect on the multiplicity ratio of the increased 
(reduced) string tension K h for kaons (protons). Thus the modification of the 
string tension by accounting for hadronic size improves the agreement between 
the computed and measured multiplicity ratios both for K~ and p. In addi- 
tion, the published results [17] for K~ without any of the above described 
model modifications are shown as dotted line. Including in the model a fla- 
vor dependent formation length, the realistic string constant and the reduced 
prehadron cross section improves our previous computation [17]. The proton 
spectra, however, cannot be described by the presented model adequately. 
We'll discuss this issue in Sec. 4.1. 

Diffractively produced vector mesons have typically large z and may influence 
the multiplicity spectra at intermediate z by rescattering prior to their decay, 
cf . Ref . [11]. A recent experimental analysis in Ref . [2] shows that diffractively 
produced vector mesons affect the pion multiplicities by a few percent only 
and therefore are neglected in our computation for all meson species. 



3 Mass-number dependence 



The dependence of hadron attenuation on the mass number A of the target 
nucleus is commonly believed to clearly distinguish the absorption and energy 
loss mechanisms for hadron attenuation. Indeed, the common expectation is 
that attenuation in absorption models is proportional to the hadron in-medium 
path length L, leading to 1— R M oc A 1 ^. On the other hand, the average energy 
loss for a parton traversing a QCD medium is proportional to L 2 , which leads 
to 1 — Rm oc A 2 / 3 . In the following we will discuss in detail the A-dependence 
for the absorption model with the main goal to obtain insights via analytic 
formulae. A numerical analysis of the A-dependence is presented in Section 4. 

For illustrative purposes and to obtain analytic formulae we take the model 
described Sec. 2 in the one step approximation where we neglect final hadron 
production in the calculation of the survival probability, since at z > 0.2 the 
average hadron formation length is longer than 10 fm, as shown in Table 1. 
To further simplify, we consider the case of a hard-sphere nucleus of mass 
number A and radius R = r^A 1 ^, with r = 1.12 fm [23]. We also consider 
only hadronization from rank 1 on, dominant in the HERMES kinematics 
except for K~ or p production. Neglecting absorption in deuterium and the 
small rescaling correction we find that the hadron multiplicity ratio Rm equals 
the hadron survival probability Sa- Therefore, the hadron attenuation 1 — R M 
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can be approximated as: 



R M PS 1 — Sa 

yrffe 2 y dyj dx——e * , (20) 



a y y tf y (/*) 

o -R(b) v 

where we omitted all flavor and hadron species dependence for simplicity. In 

Eq. (20) po is the nuclear density, r = (i^) 1/3 and R(b) = VR 2 - b 2 . With 
the mean free path of the prehadron A* = l/(cr*p ) the expression for the 
attenuation simplifies: 



l-R 



2R/{1") t r 

1 _ e ^-r) 



M = l^j(l*) 3 J dtt JdrJ due~ u 



(21) 



The attenuation is controlled by two ratios; firstly, a = (Z*)/A*, the ratio of 
formation length to mean free path and, secondly, b = 2R/(l*), the size of 
the nucleus relative to the formation length. If both quantities are small one 
can expand the integrand in powers of u, obtaining after all the integrations 
a series which starts with b 2 oc A 2 / 3 , contrary to common expectations which 
dictate A 1/3 : 

1-R M = ^ab 2 - 1(1 + a)ab 3 + 1^(1 + a + a 2 )ab* + 0[b 5 ] . (22) 

However, since a and b are of order unity for medium to heavy nuclei, the 
series converges very slowly, we cannot state whether or not it is possible to 
approximate the full expression by the leading term oc A 2 / 3 . To improve the 
convergence we fit the r-dependence of the integrand to a numerical function 
with a simpler form: 



r 

J , - , -I - < - 'M ! - ( I . _ • 



due' 

o 



£ a(u-r) 



a (I 



e- war . (23) 



l-a 

A good fit can be achieved with a constant 

w = 0.19 (24) 

independently of r and a = in the range 0<r<5, < a < 3.5 where 
the function differs substantially from zero. Using Eq. (23) in (21), one finds 
a rapidly converging expansion of the attenuation, since now the small w- 
parameter enhances the convergence: 

l-R M = \wab 2 - l(wa6 2 ) 2 + 0[(wab 2 ) 3 ] 
O i u 

= c x A 2 ' 3 + c 2 A^ 3 + 0[A 2 ] . (25) 
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Table 1 

Average hadron formation time and values of the c\ and C2 coefficients in Eq. (25) 
for 7r + production at different z values. The mass number A at which A A / 3 terms 
in (25) become comparable with A 2 ! 3 terms are also given. For each value of z, we 
have taken the appropriate average value of v measured at [2]. The large value of 
(l h ) > 10 fm justifies neglecting hadron interactions with the nucleus. 



We have computed the coefficients C\ and c 2 in Eq. (25) for 7r + -production and 
different values of z. The results are presented in Table 1. We observe that the 
series converges quickly over the whole z range. Therefore, it makes sense to 
approximate the nuclear attenuation with a power law, 

1 - R M = cA a , (26) 

where 



a 



2/3 



(27) 
(28) 



since higher-order terms will not give large corrections except at large z > 0.8. 
Indeed from Eq.(25), we see that higher order terms are negligible if 

3/2 

= A . (29) 

Values of A are given in Table 1, and become comparable to mass numbers of 
medium-heavy nuclei around z = 0.8. We will further discuss the power-law 
fit (26) in Section 4. 

We conclude with three remarks. First, a lot of information on absorption 
dynamics is contained in the proportionality coefficient c of Eq.(26), and its 
strong dependence on z needs to be taken into account when analyzing ex- 
perimental data and theoretical models. Second, the approximate power law 
1 — R M oc A 2 / 3 is not a peculiarity of the absorption model we developed, 
but has a more general validity. Indeed, it is possible to show that it appears 
whenever the probability distribution of the prehadron formation length Ay 
is analytic and has a finite limit as The origin of the A 2 / 3 is geo- 

metric and related to the integration over the 7*-quark interaction points. 
Non-analytic probability distributions may give leading powers different from 
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2/3 (cf. appendix A). Clearly, the calculation based on a hard sphere nucleus 
is not adequate for light nuclei. In the following numerical calculations of the 
full model discussed in Sect. 2 we use appropriate density profiles for all nuclei. 



4 Numerical results 

4-1 Multiplicity ratio for n, K, p and p 

The recent HERMES data have shown significant differences in the multiplic- 
ity ratios of various hadrons. As suggested in Refs. [13, 14, 16] the observed 
flavor-dependence of the multiplicity ratio may be useful to disentangle dif- 
ferent theoretical frameworks describing the hadronization process. However, 
the aforementioned models provide computations for nuclear attenuation of 
charged hadrons [14], charged pions and kaons [8, 13, 16] only. Theoretical com- 
putations for protons and antiprotons, that are expected to be very sensitive 
to the formation length mechanism as shown in Fig. 3, are given in Ref. [10]. 
These predictions generally underestimate baryon multiplicity ratio, especially 
for the krypton nucleus and for z < 0.4. 

In order to have charge-separated fragmentation functions, we use the leading 
order Kretzer parametrization [24] of the fragmentation functions for pions and 
kaons, while for protons and antiprotons we use the Kniehl- Kramer- Potter [25] 
parametrization. Predictions of our model for 7r ± , K ± , p and p multiplicity 
ratios on He, Ne, Kr and Xe are shown for the z (Fig. 5) and v (Fig. 6) 
dependence, and compared to published HERMES data for identified hadrons 
on Kr [2] and to preliminary data on He [3] and Ne [3]. (We do not include 
computations and final data on N as they are very close to Ne.) 

The z dependence reported in Fig. 5 shows nice agreement with the data for 
negative and positive pions for all the nuclei. The predictions for K + and K~ 
are in qualitative agreement with the trend shown by the data, except for the 
low-z region. This region contains contributions from both leading hadrons 
decelerated in nuclear rescattering and from secondary K + and K~ produced 
from initial pions and p-mesons that are more abundant than strange par- 
ticles. Both these contributions are not accounted in the purely absorption 
treatment of the final state interaction included in our model. The experimen- 
tal z dependence of the antiprotons data is qualitatively reproduced within 
statistical uncertainties by the model, while the proton multiplicity ratio is un- 
derestimated especially in the low z region. The discrepancy at low-z repeats 
itself in the ^-dependence shown in Fig. 6, because an average ^=0.3-0.4 enters 
in the //-dependence. A simple explanation in terms of a smaller prebaryonic 
cross-section or a larger baryon formation is clearly insufficient to account 
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for the large rise of the proton multiplicity ratio RF M above 1 at z < 0.4: 
indeed nuclear absorption can only reduce Rm below 1 and rescaling effects, 
which can increase Rm above 1 at low enough z, are too weak to describe the 
observed data. Another possible explanation for the difference in the proton 
sector is that theoretical computations do not account for the z degradation 
due to (pre)hadron rescattering nor final state interactions and decays. How- 
ever, even a full transport model like Ref. [10, 11] which takes these effects 
fully into account, fails to reproduce the rapid growth of the proton data at 
low z. In addition, the contribution of protons from target fragmentation is 
ruled out as an explanation of the effect, because the experimental momentum 
cut (\pn\ > 4 GeV) removes most of the protons knocked out of the nucleus. 
The discrepancy in the \ow-z region may also point to a non negligible diquark 
contribution as suggested in Ref. [11]. 

4-2 Mass-number dependence: a new analysis 

The simplest way of analyzing the A-dependence of hadron attenuation 1—R M 
is in terms of a power law: 

1 - R h M (z; A) = c h (z; A) A a ^ z ' A) . (30) 

For simplicity, we consider only hadron attenuation in fixed z bins, integrated 
over v and Q 2 ; a similar discussion applies for arbitrary bins in z, v and Q 2 . 
The coefficients c and a both depend in general on z and A. Indeed, as we 
have seen in section 3, the absorption model gives a power series: the lowest 
order term is A 2 / 3 , but higher order terms become important for large A and 
z. Moreover, since we expect a strong dependence of the coefficient c on z 
a single constant c cannot fit all experimental bins. Therefore, we propose 
to analyze the A-dependence of our model and experimental data alike, by 
fitting the hadron attenuation to the power law (30) leaving both a and c 
as free parameters. These parameters are strongly correlated: indeed a small 
increase of the exponent a can be compensated by a small decrease of the 
coefficient c, and vice-versa. Therefore the fit results are best presented as x 2 
contour plots in the (c, a)-plane, showing the aforementioned correlations. 

The best-fit parameters c and a are determined by chi-square minimization, 
i.e. the x 2 merit function 

X 2 (c, a)=J2 *[(1 - R h M )(A t ) - cAf] 2 (31) 

i i 

is minimized with respect to c and a. Here <7j is the uncertainty of the theo- 
retical or the experimental points respectively. For the experimental data we 
use the statistical uncertainty only, while in the case of the theoretical com- 
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putations we use the quadratic sum of the uncertainties corresponding to the 
precision of the numerical computation («1%), and to the choice of model 
parameters. In our model we fit the value of the prehadron-nucleon cross sec- 
tion u ^ to experimental data for n + production on Kr. Hence the theoretical 
relative uncertainty is determined by the experimental relative uncertainty 
of 7r + attenuation on Kr, which approximately equals 6% independent of z. 
Then, we assume that the uncertainty of <t* similarly yields a 6% relative un- 
certainty on hadron attenuation independently of z and A, as well. A change 
in the other parameters of the model would be reflected in a slightly different 
value of o"*, but would not affect the estimate of the theoretical uncertainty of 
hadron attenuation. 

The x 2 contour plots in the (c, a)-space are computed as constant x 2 bound- 
aries, enclosing a region such that 

Xmin *^ X ^ Xmin + A X 2 (32) 

where Xmin * s the minimum of the x 2 function, obtained at the best-fit parame- 
ters. We fix Ax 2 = 4.61 which corresponds, in the case of normally distributed 
fit parameters, to the region which covers 90% of the total probability distri- 
bution. 

The dotted x 2 contours in Fig. 7 represent the fit to the numerical computation 
in various z-bins of the full theoretical model as it is described in section 2, 
including 4 He, ( 14 N), 20 Ne and 84 Kr nuclei. In order to have a direct comparison 
with the theoretical mass number dependence analysis of the pure absorption 
model in section 3, the rescaling of parton distribution and fragmentation 
functions is not included in the numerical computation yielding the dashed 
contours. The solid x 2 contours are obtained from the published HERMES 
data on ( 14 N [1]) and 84 Kr [2] and the preliminary data on 4 He [3] and 20 Ne [3]. 
The notation ( 14 N) indicates that this nucleus is included in the fits of theory 
and experiment only for z > 0.55, because there are no experimental data in 
lower z-bins. The fits at z < 0.45 contain only three data points, namely He, 
Ne and Kr, and 2 fit parameters. Therefore, the fit yields in these bins at best 
a rough estimate only. Hence, we consider the results in these z bins to be less 
reliable. 

The presented results on He, (N), Ne and Kr at z > 0.55 for the pure absorp- 
tion model computation are qualitatively in agreement with the trend shown 
by the experimental data, except for the z = 0.85 bin, where the experimental 
fit gives the largest a due to 1 — R M < value for 4 He. At small z, the full 
model contours behave similarly to the pure absorption contours, but the two 
become more and more separated as z is increased. Furthermore, the full model 
increasingly disagrees with experimental data when z > 0.75. This shows the 
power of the proposed kind of analysis, which is a possible tool to disentangle 
different theoretical models. 
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Table 2 

Centroids of the contour plots in Figs. 7 and 8 with their uncertainties for the fit 
1 — Rm = c(u, z, h)A a at fixed z bins, both for the experimental data and the pure 
absorption model calculation. The nuclei included in the fits are shown in the table, 
the parenthesis on (N) indicate that nitrogen is included in the fit only for z > 0.55. 
Note, that the theoretical results include different sets of nuclei. 



The slope of the \ 2 contour plots shows a trend to decrease by increasing z. In 
general, the slope of the contours in the (c, a)-plane for the different z-bins is 
estimated by calculating the variation of 1 — I& M = c(A) a at an average mass 
number (A): 

A(l- R h M ) = {A) a Ac + ca{A) a ' 1 Aa. (33) 
A change Ac of c is correlated with a corresponding change Aa of a given by 

Aa = -^-Ac . (34) 
ca 

This equation gives the negative slopes of the contours. As the product of c 
and a generally increases with increasing z, this implies a decreasing slope. 

The best-fit parameters, with their uncertainties are given in Table 2, both for 
the fits of the experimental data and the pure absorption model predictions. 
By showing the best-fit parameters of the pure absorption model computation 
we can quantitatively compare deviations from the leading order A 2 ^ 3 depen- 
dence predicted in section 3. The uncertainties of the best-fit parameters are 
determined by projecting the contour plots on the c- and cc-axes respectively. 
Therefore, here as well, the parameter uncertainties correspond to a joined 
90% confidence interval in the case of normally distributed parameters. 
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The values of a derived from both the pure absorption model fits and the 
experimental fits on He,( 14 N), Ne, Kr at z > 0.55 are compatible, within their 
uncertainties, with A 2 / 3 behavior. Deviation from this behavior should be 
observed for heavier nuclei than krypton and at high z as shown in Table 1. 
In order to observe the breaking of the power law predicted by the model 
and to increase the statistical significance of the lowest z-bins, we perform 
the same computation by including xenon nucleus. Experimental data on Xe 
have been recently collected by the HERMES experiment, but the analysis 
is as yet in progress. The theoretical contour plots are shown in Fig. 8 and 
the best-fit parameters for the pure absorption model with their uncertainties 
are summarized in Table 2. Since we include a nucleus with a mass number 
comparable to A for large z bins, as shown in Table 1, the power law A 2 / 3 is 
broken significantly producing a reduction of the a values. 

In Fig. 9 the published HERMES data for R M for 14 N [1] and 84 Kr [2] and 
the preliminary data for Rm for 4 He [3] and 20 Ne [3] are displayed as 1 — Rm 
for two selected z-bins (z = 0.45, 0.75) with full diamonds and triangles. The 
empty symbols show the corresponding results of the pure absorption model 
for the above nuclei plus 131 Xe. In addition we plot with solid (dashed) lines 
the best fits to 1 — R M = cA a of the theory results including He,N,Ne,Kr 
(He,N,Ne,Kr,Xe). One observes that the inclusion of xenon flattens the curve 
for the attenuation 1 — Rm for large values of z. The exponent extracted by the 
proposed (c, a) fit can be regarded as an average exponent on the considered 
interval of atomic masses. Beyond A m 80 and for large values of z the xenon 
point lies below the A a curve fitted to all nuclei. 

In order to investigate a breaking of the A 2 ! 3 power law at large A we show 
Fig. 10, where the z-dependence of theoretical a values in our pure absorption 
model is compared with the z-dependence of a values derived from experi- 
mental data. The xenon nucleus shifts the upper band limit of the theoretical 
predictions to lower a values. But the a values with Xe are compatible with 
the a values without the Xe nucleus over the whole z range. Nevertheless the 
A-dependence of hadron attenuation in nuclei is a promising tool in order to 
distinguish different theoretical assumptions. 



5 Summary and conclusions 

Rescaling of the structure and fragmentation functions together with (pre) 
hadronic absorption describe both HERMES and EMC data on the nuclear 
modification of hadron production in DIS processes [17]. This paper improves 
the original model of Ref. [17], by correctly evolving quarks, prehadrons and 
hadrons without factorizing the absorption process and hadronization pro- 
cesses. The hadron formation length is shown to be flavor dependent. Espe- 
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cially, negative kaons (K~ = su) as well as antiprotons (p = uud) cannot 
be formed by a struck valence quark which is the dominant contribution at 
HERMES energy. This yields a smaller formation length for negative kaons 
and antiprotons implying more hadron suppression compared to positive kaons 
or protons respectively. Furthermore the so far neglected transverse extension 
of the produced hadrons leads to slightly different formation lengths for dif- 
ferent hadron species. We now use a reasonable value for the string tension 
k — 1 GeV/fm, and extract the prehadronic cross-section a* = 2/3cr h from 
data on pion production on krypton. This shows that the prehadron does not 
have the full hadronic size yet. 

The presented model correctly describes pion and kaon multiplicity ratios. A 
different magnitude of the multiplicity ratios for positive and negative kaons 
originating from a flavor dependent formation length is reflected in the experi- 
mental data. In the baryon sector, although computed antiproton mutliplicities 
agree with experimental data inside their uncertainties, model computations 
disagree with proton data especially at small z. This discrepancy may be as- 
cribed to a non negligible diquark fragmentation contribution in that region. 
Prehadron rescatterings and final state interactions are insufficient to fully 
account for the effect. 

We have applied our model to the analysis of the mass number dependence 
of the hadron attenuation 1 — Rm- Contrary to the common expectation we 
found a cA 2 ^ 3 behavior in leading order, with higher order corrections increas- 
ing in magnitude with increasing z and A. Hence, one must not expect a strict 
power law, but rather an effective behavior cA a . Furthermore, the pro- 
portionality coefficient c depends strongly on the kinematic variables, which 
rules out a simple analysis with a globally fixed c. Therefore we have proposed 
to analyze the A-dependence of hadron attenuation in terms of a power-law 
fit 

l-R M = c A a (35) 

with both c and a as fit parameters. Their correlations, displayed as x 2 con - 
tours in the (c, a)-plane, are quite sensitive to the underlying model assump- 
tions and may disentangle different theoretical models. 

Qualitative agreement is found between the theoretically and experimentally 
determined x 2 contours, with a ~ 2/3. Contrary to common expectations, 
a pa 2/3 also in the presented absorption model. This value is not a peculiar- 
ity of our model, but a common feature of a quite large class of absorption 
models (a sufficient requirement being an analytic probability distribution in 
the formation length Ay, with a finite limit as Ay— >0). This fact bars the much 
advertised use of the exponent a = 2/3 as a direct experimental evidence of 
parton energy loss following the QCD "L 2 " law. 
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In different models the breaking of the a = 2/3 power law happens at different 
values of A and z. Therefore collecting data up to heavy nuclei like Sn and Pb, 
which is in progress at Jefferson Laboratory with the CLAS experiment [26], 
will help unraveling the space-time dynamics of the hadronization process. 
The proposed analysis of the A-dependence of the hadron attenuation may 
thus provide a new method to differentiate between absorption and energy 
loss effects. In this respect we encourage the authors of competing models to 
perform a similar (c, a) analysis to have a full set of theoretical results when 
the new data from HERMES and CLAS measurements will be available. 
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A Generic absorption model 



In this appendix we show that hadron attenuation behaves as 1 — Rm = 
cqA 2 / 3 + 0(A 1 ) for a large class of absorption models and is not a peculiarity 
of our model. A sufficient requirement is that the probability distribution of 
the prehadron formation length Ay is analytic and has a finite limit as Ay^O. 

We consider the generalized 1 step approximation described in Sect. 2 by 
substituting in Eq. (20) the exponential probability distribution for prehadron 
formation with a generic distribution V: 

-e W _> p _ JL) , (A.l) 



(/*) A V A 

where A is a typical scale of the fragmentation process, and V(u) is normalized 
to 1: 



/ duV{u) 
Jo 



1 . (A.2) 
The absorption factor (21) is then generalized as follows: 

2R/A t r 

1-Rm = ^A 3 J dttjdrj duV{u) [l - e£ (M ~ r) ] (A.3) 



(note that integration over t and u is the integration over all possible 7*g 
interaction points). This generalization captures the essential features of the 
A-dependence of the hadron suppression factor 1 — R M in most absorption 
models found in literature. A few examples will be discussed below. 

A quite general probability distribution can be defined as follows. Let's V 
be an analytic function on (0, M) with M > or M = oo. Let's n be the 
smallest real number such that lim u ^ u~ n V(u) = p with p ^ 0, oo. Then 
the function 



/(«) 



u- n V(u) u G (0,M) 
Po u = 



can be expanded in powers of u around u — 0: 

f(u) =p + Pl u + O(u 2 ) . (A.5) 
This means that 

V{u) = Po u n + 0{u n+l ) . (A.6) 
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By the normalizability requirement (A. 2) the exponent n must satisfy 



n > -1 . (A.7) 

However, to our knowledge, there is no probability distributions considered in 
literature with n < 0. Now, if we expand the integrand of (A. 3) in a Laurent 
series around u — and perform all the integrations, we obtain 

1 D vrpo p 2 n+5 A 4 iT+ 5 t ^[R n+6 \ 

l — n M = -r-r-/ r~7T7 TT^T 1 ^T7 ~^T~T~ » s + U 



2A (n + l)(n + 2)(n + 3)(n + 5) A* A n + 5 \A n + 6 



(A.8) 



If we now use R = r^A 1 ^ 3 we obtain 

1-R M = c Af +§ + 0(A* +1 ) , (A.9) 

where the constant Co is easily read off Eq. (A.8). We see that the exponent of 
the first-order term is always larger than 1/3 by virtue of Eq. (A.7). However, 
since no models in literature have n < 0, the exponent is in practice always 
larger than or equal to 2/3. From Eq. (A.9) we see that in absorption models 
such that lim u ^ 'P( u ) — Po > 0, which have n — 0, the expansion of 1 — Rm in 
powers of A starts with an A 2 / 3 term, as claimed. The exponential distribution 
considered in this paper falls in this class of models. 

An example of a model with n > is the bremsstrahlung model of Ref. [13]. In 
that paper the computed probability distribution satisfies V(u)^0 as u— >0, 
so that n > and 

(1 - Rm)[i 3] = c [m A a w + Q(A a ^) (A.10) 

with ari3] > 2/3. However, except for very small u or small z, the probability 
distribution is well approximated by an exponential distribution. Hence, for 
realistic finite-size nuclei we expect higher order corrections to be important 
and to make the A dependence of 1 — R M effectively close to A 2 / 3 . 

Models with non-analytic probability distribution can have different behaviour. 
For example, in Ref. [11] all prehadrons are assumed to be formed immediately 
after the interaction of the virtual photon with a quark (though prehadrons 
created from the middle of the color string are assumed not to interact with 
the medium). Therefore, absorption of rank 1 prehadrons can be modeled as 
in Eq. (A.3) with 

V(u) = 5(u) . (A.ll) 
Carrying out all the integrations in (A.3) one finds 

(l- J RAf)[ii] = c [ ii ] ^ 1 / 3 (A.12) 
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with no higher order corrections. However, this A 1 / 3 power law may change 
when also considering the other processes included in the model, like diffractive 
meson production (important at large z) and higher-rank prehadron produc- 
tion, elastic (pre)hadron rescatterings, resonance decays (important at smaller 

z). ' 

On the other hand, a non-analytic probability distribution is not a guarantee 
for an A 1 / 3 power-law. Indeed, the Bialas-Chmaj model (BC) of Ref. [6] can 
be shown [17] to reduce to Eq. (A.3) with A = (/*) and 

V{u) = 5[u - (1 - e- r )] . (A.13) 

Carrying out all the integrations in (A.3) one finds 

(l-R M ) m =c [6] A 2 / 3 + o(A 2 / 3 ) . (A.14) 

This concludes our proof that an A 2 / 3 power law for the hadron suppression 
factor of not too heavy nuclei is a general feature of a quite large class of 
absorption models. 
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Fig. 5. Computed multiplicity ratios (cf. Eq. 1) for pions, kaons, protons and an- 
tiprotons as a function of z for He (dashed line), Ne (long-dashed line), Kr (solid 
line) and Xe (dot-dashed line). The HERMES data on Kr (solid circles) [2] and the 
preliminary data on He (closed triangles) and on Ne (open circles) [3] are shown 
with their statistical uncertainties. To improve the readability of the figure, He and 
Ne experimental points have been shifted by 0.015 to the left and to the right, 
respectively. 
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Fig. 6. Computed multiplicity ratios (cf. Eq. 1) for pions, kaons, protons and an- 
tiprotons as a function of v for He (dashed line), Ne (long-dashed line), Kr (solid 
line) and Xe (dot-dashed line). The HERMES data on Kr (solid circles) [2] and the 
preliminary data on He (closed triangles) and on Ne (open circles) [3] are shown 
with their statistical uncertainties. To improve the readability of the figure, He and 
Ne experimental points have been shifted by 0.25 GeV to the left and to the right, 
respectively. 
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Fig. 7. x 2 contours for the fit 1 - R M = c{v, z, h) A a on 4 He, ( 14 N), 20 Ne and 84 Kr 
nuclei, in fixed z-bins. The fits to our model computation are shown by the dashed 
X 2 contours in the case of the pure absorption model and by dotted \ 2 contours in 
the case of the full model (absorption + rescaling). The solid contours show the fits 
to the HERMES pion multiplicity ratios on (N [1]), Kr [2] and to the preliminary 
data on He and Ne [3] , accounting for statistical uncertainties only. The parenthesis 
on (N) indicate that nitrogen is included in the fit only for z > 0.55. Note that in 
the last z bin the contour for the full theory computation is partially out of the plot 
range. 
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Fig. 8. x 2 contour plots for the fit 1 — Rm = c(v, z,h)A a on the pure absorption 
model (dashed) and the full model (dotted) computation for 4 He, ( 14 N), 20 Ne, 84 Kr 
and 131 Xe nuclei, in fixed z-bins. Note that the contour for the full theory compu- 
tation in the last z bin is outside of the plot range, to the right. 
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Fig. 9. The published HERMES data for R M for 14 N [1] and 84 Kr [2] and the 
preliminary data for Rm for 4 He [3] and 20 Ne [3] are shown as 1 — Rm for z = 0.45 
and z = 0.75 as filled diamonds and filled triangles respectively. The pure absorption 
model results for 1 — Rm in the same z-bins and for the same nuclei plus 131 Xe are 
shown by empty symbols. Note that we have shifted the absorption model values 
slightly to the right to avoid overlap with the experimental points. The solid lines are 
a best fit of the pure absorption model results to 1 — Rm = cA a including values for 
A = 4, 14, 20, 84 (He,N,Ne,Kr) and the dashed lines represent a fit including values 
for A = 4, 14, 20, 84, 131 (He,N,Ne,Kr,Xe). 



29 



i — i — i — i — i — i — i — i — i — i — i — ] — i — i — r 




Q I i i i I i i i I i i i I i i i I 

0.2 0.4 0.6 0.8 1 

z 

Fig. 10. Values of a as a function of z derived from experimental data (dots) and 
our pure absorption model computation (bands). The nuclei included in the fit 
are shown in the legend (the HERMES data on He and Ne are preliminary). The 
notation (N) indicates that this nucleus is included in the fit only at z > 0.55. 
Note that experimental errors are uncorrelated, but theory errors are point-to-point 
correlated. 
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